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ABSTRACT 

Understanding the non-adiabatic pressure, or relative entropy, perturbation is crucial for stud- 
ies of early-universe vorticity and Cosmic Microwave Background observations. We calculate 
the evolution of the linear non-adiabatic pressure perturbation from radiation domination to 
late times, numerically solving the linear governing equations for a wide range of wavenum- 
bers. Using adiabatic initial conditions consistent with WMAP seven year data, we find nev- 
ertheless that the non-adiabatic pressure perturbation is non-zero and grows at early times, 
peaking around the epoch of matter/radiation equality and decaying in matter domination. At 
early times or large redshifts (z — 10, 000) its power spectrum peaks at a comoving wavenum- 
ber k « 0.2/i/Mpc, while at late times {z = 500) it peaks at A: w 0.02/i/Mpc. 
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1 INTRODUCTION 

Our understanding of the physics of the universe is improving 
steadily. Recent observational successes, in particular Cosmic Mi- 
crowave Background (CMB) experiments and Large Scale Struc- 
ture (LSS) surveys ( |Komatsu et al.|2011a[|Ade et al. 201 1 ; Abaza- 
[jian et al.|2009| l, strongly favour the cosmological standard model, 
in which the anisotropies of the CMB and the LSS are sourced by 
quantum fluctuations formed during a period of accelerated expan- 
sion in the early universe - inflation ( [Liddle & LythpOOO ). 

Since the decay of the scalar field(s) driving inflation into the 
standard matter fields is not fully understood, the power spectrum 
of the fluctuations in the scalar field is mapped onto the spectrum 
of the comoving curvature perturbation, TZ, at horizon crossing. 
Since TZ is conserved on large scales for adiabatic perturbations, 
the primordial power spectrum is not affected by the small scale 
physics, and can be used to set the initial conditions for Boltzmann 
codes l|Bertschinger|1995||Seljak & Zaldarriaga|1996[ [Lewis et al.| 
|2000 : Doran 2005, Lesgourgues 20fl]l. The current data allow for 
a small non-adiabatic contribution to the primordial perturbations, 
only providing an upper bound ^Komatsu et al.|201 lb| l. 

This contribution is described by the non-adiabatic pressure, 
or entropy, perturbation, SPnad- To date this has been calculated 
mainly to get a handle on the evolution of the curvature perturbation 
( [Wands et al.|2000) , and hence the primordial power spectrum, and 
also played a key role in characterising the initial conditions ("adi- 
abatic" versus "isocurvature" ( [Liddle & Lyth|2000l l). Observations 
of the CMB are fully consistent with purely adiabatic initial condi- 
tions jKomatsu et al.|201 Ib^ , although a level of primordial isocur- 
vature can be allowed (references i Valiviita & G iannanto nio|2009| 
[Mangilli et al.|2010[[iret al.|2011| i provide an non-exhaustive ex- 
ample of recent studies). 



However, after the initial conditions have been imposed, there 
are no restrictions on the non-adiabatic contribution. In this let- 
ter we focus on the relative entropy perturbation, which naturally 
arises in any multi-component system. This contribution to 5Pnad 
has so far not been studied in detail in realistic models, and we 
calculate its evolution numerically, starting from adiabatic initial 
conditions in agreement with the WMAP7 concordance model. 

A non-zero entropy perturbation sources the evolution of the 
curvature perturbation, which will not be conserved on large scales 
beyond the time when the initial conditions are set. In practice, this 
will affect any calculation assuming conservation of the curvature 
perturbation after z ~ 100, 000. Furthermore, it has been shown 
that the non-adiabatic pressure perturbation can generate vorticity 
in the early universe (Christo pherson et al.[[2009) . Although this 
vorticity is generated at second order in perturbation theory, it 
might act as foreground and contribute to the primordial B-mode 
polarisation signal in the CMB on small scales. 



2 DEFINITIONS 

We consider linear, scalar perturbations to a flat Friedmann- 
Robertson- Walker (FRW) spacetime, following [Malik & Wands [ 
( [2005) 1. The total density and the total pressure including back- 
ground and perturbations, p and P, respectively, are related to the 
density and pressure of the component fluids p = Pa , P = 
X^tj ' where Greek letters label the individual fluids. 

The total pressure perturbation can be split into an adiabatic 
and non-adiabatic part 5P = 5Pnad + ci5p , where ci is the adia- 
batic speed of sound, Cg = P / p , and an overdot denotes a deriva- 
tive with respect to conformal time, 77. 
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In the presence of more than one fluid, the total non-adiabatic 
pressure perturbation, 5Pnad, may be further decomposed |Kodama| 
|& Sasaki! 1984) 1 as JPnad = (5-Pintr + 5-Prci • The first part is the 
sum of the intrinsic entropy perturbation of each fluid, which van- 
ishes for barotropic fluids. We assume here for simplicity that all 
individual fluids have a constant equation of state and are hence 
barotropic. 

The second part of the non-adiabatic pressure perturbation is 
due to the relative entropy perturbation 5^/3 between different flu- 
ids, defined as 



5c, fl 



spa SpjS 



pa 



(1) 



where T-i is the conformal Hubble parameter, T-i = a. /a. The rela- 
tive non-adiabatic pressure perturbation between multiple fluids is 
then 



rol = - ^ PaP0 (4 - 4) "^"/S 



(2) 



a,0 



P a, 13 

with = Pal pa the adiabatic sound speed of each individual 
fluid. It is straightforward to see that the non-adiabatic pressure 
perturbation is invariant under a gauge transformation. This also 
implies that if it is non-zero in one gauge, it will be non-zero in any 
other gauge; i.e. it can not be "gauged away". 

By definition, the relative non-adiabatic pressure perturbation 
does not contribute to the individual fluid evolution equations and 
instead enters the Einstein equations. It is perhaps in this context 
that the non-adiabatic pressure perturbation is most familiar, gov- 
erning the evolution of the curvature perturbation on uniform den- 
sity hypersurfaces, on large scales jWands et al.|2000y 

Boltzmann codes typically avoid the Einstein equations in 
which the non-adiabatic pressure perturbation would enter - these 
have instead been used as consistency checks - and so the system 
we are evolving remains closed at linear order. 

However, at higher orders in perturbation theory, knowledge 
of the behaviour <5Prci becomes important. In particular, the 
non-adiabatic pressure produced by linear perturbations will have a 
significant impact on the curvature perturbation on uniform density 
hypersurfaces, C,, at a non-linear level. A first study in this direction 
was undertaken at second order by |Malik| ( |2005[ l, showing the 
dependence of C,2 on the relative entropy perturbations. It has 
also been shown that the non-adiabatic pressure perturbation 
can generate vorticity in the early universe ( [Christopherson et al.| 
|2009 t. Furthermore, if a modification to a Boltzmann code includes 
a perturbed total pressure the non-adiabatic pressure perturbation 
should be included. 



3 DYNAMICS AND INITIAL CONDITIONS 

The concordance cosmology contains five fluids: baryonic matter, 
cold dark matter, photons, neutrinos and some form of dark en- 
ergy. We assume the neutrinos to be massless. Defining the equa- 
tion of state parameter Wa = Pa/ pa, baryons (6) and CDM (c) 
have Wb = Wc ~ ci ~ ci ~ Q- Photons (7) and neutrinos {v) 



are relativistic species and have w-, 



1/3. 



Photons and baryons are coupled together by Compton scattering, 
while neutrinos stream freely. We employ a flat ACDM cosmology 



consistent with the WMAP 7-year results ( [Komatsu et al.|2011b| l 
with the parameters 



fit/i^ = 2.253 X 10'^ , n^h'^ = 0.112 , 
= 0.728, ft = 0.704. 



(3) 



Although we employ a cosmological constant, our approach also 
allows for a dark energy with 7^ 1 and therefore dark energy 
density perturbations. These would contribute to the non-adiabatic 
pressure perturbation at late times. Since our focus is on the early 
universe where dark energy perturbations are expected to contribute 
negligibly, results from ACDM will be very similar to those from 
more general models. Of course, formally we should include dark 
energy perturbations iPark et al. 2009 Christopherson 20101, and 
these will be investigated in future work. 

Since the non-adiabatic pressure perturbation is gauge- 
invariant it can be evaluated in any gauge, and we choose for 
numerical simplicity the synchronous gauge, in which the line- 
element takes the form 



ds 



(4) 



employing the conformal time r]. The residual gauge freedom 
is removed by choosing to co-move with cold dark matter, 
i.e. Vc{k,r)) — 0. The go verning equations are then presented in 
Ma & BertschingerHl995] n 



Taking the entropy perturbation between each species to van- 
ish at some early epoch gives the fluids the initial conditions in 
tight-coupling and deep radiation domination \Ma & Bertschinger] 
[T995I 1 



S-y — 5,y = ^Sb = |(5c ~ 



3 



-V. 



1 ■ 



0. 



(5) 



These are supplemented by initial conditions for ?/>, E and the neu- 
trino shear which we do not present here; for clarity we have also 
neglected small corrections due to the non-zero matter density even 
in the early times, although they are included in the code. These ini- 
tial conditions are implemented at Zi ~ 100, 000 to ensure that all 
linear modes are strongly superhorizon when the code begins. The 
constant = pu/ {p~, + p^) is the relative neutrino abundance. 

The system of equations can now be integrated to find the fluid 
and metric perturbations at arbitrary time. The power spectrum of 
a perturbation (5(k) is given by 



Ps{k,i^) = <<5(k)r(k')) 



Pc(fc)l<5(fe)l'(27r)='5^(k-k') 



(6) 



where Vc{k) = yl*(fc/fc*)"="^ is defined using the curvature per- 
turbation ^ TZ which is constant on superhorizon scales. The 
amplitude A* at the pivot wavenumber fc* take the WMAP7 values 
A* = 2.42 X 10"'-* at it* = 0.002/iMpc"\ Figure [I] shows the 
power spectrum of the baryon density contrast Sb{k), 



2n^ 

Pb{k,v)^-^Vdk)Mk)\ 



(7) 



for the range of wavenumbers and redshifts we sample. We note 
that the maximum plotted redshift Zmax ~ 50, 000 and is therefore 
significantly later than the initial redshift Zi ~ 100, 000. 



^ In terms of the variables h and rjMB in|Ma & BertschingerHl995|, the 
spatial curvature perturbations ^ and E (in Malik & Wands ( 2005 1) are 
■0 = r^MB and_E = — (/i + 6»)mb)/(2A;^), while the velocity perturbation 

is Va = Oa/k'^- 
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Figure 1. Power spectrum of the baryon density contrast, Pb{k, rf). 



4 RESULTS 

We can now derive an analytical approximation for the non-adiabtic 
pressure perturbation. Defining the scaled entropy differences be- 
tween two fluids fluids by 



Aq/3 = (1 + VJlj)5a — (1 + 'Wa)Sfi , 

and using the fact that the density contrasts evolve as 



5c + \h ■- 



5^ + |/i + \k^v^ 



0, Sb + ^h + k^vb^Q, 



0, 



3 ' 3 

we find that the scaled entropy differences are governed by 

A^ft = Ifc^ {vh — v^) , A^c = 



A^6 



3 



A,c 



" 2 A/ . 



(8) 



(9) 



(10) 



These equations are gauge-invariant, and apply at all times. 

From Eq. l|5j we see that on superhorizon scales, in radia- 
tion domination and during tight-coupling, the CDM terms will 
dominate. Furthermore, the present neutrino abundance implies 
Rv ~ 0.4 and so 



A„c > A-,c > A^i, > A 



-fb • 



(11) 



Assuming tight-coupling to be exact, implying A^t ~ 0, integrat- 
ing Eqs. \10\ in the early universe then gives 

3/216 xHgn.(15 + 12i^.) ^,Ya--^fc^a. (12) 
BttGI 1+ |(15 + 4i?,) 



In finding this we have used that in radiation domination a oc 77. A 
rough approximation to the behaviour in early matter domination 
can be found by setting a oc rf , giving the rapidly decaying 

5P,^iOik'^a''^{l + mMa/AQ.Ry'^ (13) 

which suggests that SP^ci reaches a maximum at or around the 
epoch of matter/radiation equality. However, Eq. \\2i) applies only 
very close to equality as it assumes a similar growth of pertur- 
bations to that in radiation equality; the most we can say is that 
we expect the non-adiabatic pressure perturbation to decay early 
in matter domination. This will be confirmed by a more accurate, 
numerical analysis. 

Given the primordial power spectrum ^{(fc) oc k"^^^ the 
power spectrum of the non-adiabatic pressure perturbation in the 
early universe is then 
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Figure 2. Power spectrum of the non-adiabatic pressure perturbation 

-P«P„l{ft,77)- 



The non-adiabatic pressure perturbation on large scales and at early 
times grows with scale factor. It is also strongly dependent on k and 
decreases rapidly on superhorizon scales, which is to be expected. 

To calculate the evolution of the non-adiabatic pressure per- 
turbation we employ a modified version of the CMBFast code jSel-| 
^ak c& Zaldarriaga 1996 ; Bertschinger 1995 1 to recover the velocity 
perturbations for the full system from radiation domination to the 
present day, across a wide range of scales. The non-adiabatic pres- 
sure perturbation is then found by integrating Eqs. \\Q\ . 

The power spectrum of the non-adiabatic pressure perturba- 
tion is plotted across scale and redshift in Figure |2] The oscilla- 
tions before recombination are prominent, and SPici indeed peaks 
at approximately matter/radiation equality. It is orders of magni- 
tude more significant on smaller scales than on larger scales and 
at early times peaks at approximately k ~ 0.1/i/Mpc. Across all 
wavenumbers SP^^i decays rapidly after matter/radiation equality. 
The non-adiabatic pressure perturbation is negative. 

The left panel of Figure |3] shows the evolution of this non- 
adiabatic pressure perturbation for wavelengths between k — 
lO^^ftMpc"^ and k = 10~^/iMpc. The growth with scale-factor 
at the earliest times is apparent, as is the decay during matter dom- 
ination. However, the period at which SP^d peaks is extended, be- 
ing approximately within an order-of-magnitude of the peak be- 
tween z ~ 10, 000, significantly before matter/radiation equality, 
and recombination at z ~ 1, 100, after which it decays rapidly. 
For very short wavelengths SP^ci peaks at earlier times and decays 
more rapidly. The decay after recombination goes as Psp^^^ oc 
for smaller scales and as Psp,^^ oc for larger scales, imply- 
ing SPrci oc a~'^ and SPrc\ oc a~^^^ on small and large scales 
respectively. 

The right panel of Figure [3] shows the dependence of the non- 
adiabatic pressure perturbation on wavenumber for a range of red- 
shifts from deep radiation domination to the present epoch. The 
decay towards the present epoch is clear, and the scale at which the 
perturbation has most power grows from k ~ 0.2 in deep radiation 
domination to ~ 0.02 at a redshift of 2 « 500. The power spec- 
trum scales as Psp^^^i oc k~^^'^ on small scales and as Psp^^i oc k^ 
on large scales. The non-adiabatic pressure perturbation itself then 
scales as SP^ci oc fc^''* on small scales and SP^ei oc k^ on large 
scales. 

These results can be summarised as 



(5P.ei(k)5P*i(k')) oc k^ 



(14) 



SPial OC a 



k*, fc<10"^/i/Mpc 
fci/^, fc>10"^/i/Mpc 



(15) 
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Figure 3. Power spectrum of the non-adiabatic pressure perturbation Psp^^^{k,r)) ^s a function of redshift for set wavenumber (left); and as a function of 
wavenumber for set redshift (right). 



at early times, and 



5P, 



a~ 



k < 10"^ - 10-*;i/Mpc 
k > 10"^VMpc 



(16) 



at late times. This contrasts with |Christopherson et al.| pOl 
which assumed 5Pnad oc k"a~^ with a 7^ 1 before setting a = 2. 
This choice underestimates the gradient and overestimates the rate 
of decay of 5Pnad on large scales while failing to model the be- 
haviour on small scales. We would like to emphasise that we are 
working with small but finite wavenumbers; in the limit fc — > the 
non-adiabaticity vanishes, and the rapid decay (~ k*) on super- 
horizon scales is very clear. 

The contribution to the non-adiabatic pressure perturbation 
due to the entropy difference between fluid species a and /3 is 



rcl,Q:/3 — 



PaP/3 ((1 + 1^/3)60, - (1 + Wa)5fj) 

3{P + P) 



(17) 



Figure |4] shows the contributions of these terms at wavenum- 
bers k = W'^hMpc'^ and k = 10~^/iMpc~\ At early 
times, the CDM/neutrino contribution dominates slightly over 
the CDM/photon contribution, as was expected in Eq. 
The baryon/neutrino term is significantly smaller, and the 
baryon/photon term (as expected) is smaller again. 

After recombination, this picture changes. The decoupling 
of the photons and baryons causes their relative entropy to in- 
crease significantly until at 2; ~ 100, (5Prci.76 SP^cLvb- 
At late times the baryon perturbations and the CDM perturba- 
tions become approximately equal, while v^, and behave sim- 
ilarly to one-another. It is no surprise that the CDM/neutrino 
and CDM/photon contributions exhibit the same behaviour as 
the baryon/neutrino and baryon/photon terms, and so after re- 
combination, SPj-ci^jc > SPrci,iyc- The difference between the 
CDM/radiation and baryon/radiation terms is then due to the rel- 
ative abundance of CDM over baryons; in the concordance model, 
the CDM terms dominate by a factor of 4 over the baryon terms. 

On small scales the behaviour is similar to the large-scale be- 
haviour, except that 5Prci,7c -Sj SPid.vc for the entire period sam- 
pled. For 2; > 5 X 10* this would be expected to be reversed. 
For all times and across all linear scales, one can therefore find 
a reasonable approximation to the non-adiabatic density perturba- 
tion by considering only the relative entropies between the radiative 
species and the cold dark matter. 



5 DISCUSSION 

We have calculated the non-adiabatic pressure perturbation in the 
WMAP7 concordance cosmology, finding both its magnitude and 
its scaling properties. It is non- vanishing and grows in the early uni- 
verse, becoming most significant at around matter/radiation equal- 
ity. For a wide range of wavenumbers it remains within an order 
of magnitude of its peak until recombination after which it decays 
rapidly. 

Although the existence of the non-adiabatic pressure perturba- 
tion does not influence existing first-order CMB calculations em- 
ploying Boltzmann codes such as CAMB or CMBFast, its study is 
nevertheless extremely important. Firstly, this quantity is not actu- 
ally calculated in Boltzmann codes and yet has a concrete phys- 
ical meaning - it is directly related to the entropy perturbation. 
One example we used to motivate its importance is the vorticity 
it will produce at second-order, the study of which is ongoing, and 
the magnitude of which could be sufficient as to provide an ad- 
ditional foreground to the CMB polarisation signals. In turn this 
may produce magnetic fields that could influence recombination 
physics or provide seed fields for the observed cluster magnetic 
fields, or to influence large-scale structure formation in the more 
recent universe. Secondly, even though the non-adiabatic pressure 
is formed of linear combinations of quantities evaluated in Boltz- 
mann codes, it cannot be accurately recovered by naively combin- 
ing these quantities. We have assumed entirely adiabatic initial con- 
ditions. On the largest scales and at early times, then, we have for 
instance ~ (4/3)(5(,- The non-adiabatic pressure depends on the 
difference between these two quantities, which is extremely small. 
Calculating this numerically is then rather unstable; naively recon- 
structing the non-adiabatic pressure simply from the perturbations 
in one of the Boltzmann codes could give entirely the wrong be- 
haviour on very large scales and at very early times[^As such, one 
needs to modify the codes to evaluate the entropy differences di- 
rectly, as we have in this work. 

This perturbation will have various effects. Whereas the cur- 
vature perturbation is conserved on large scales at early times while 
the universe can be modelled by a single barotropic fluid, this is no 
longer case at later times when the universe has to be modelled as 
a multi-fluid system. Another effect is the generation of vorticity 



^ We found an extreme case where the reconstructed non-adiabatic pres- 
sure decays at early times, since the rounding error in 5-y — (4/3)<5i, grows 
ever more dominant for increasing redshifts. 
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at second order. The assumptions in [Christopherson et al.|pbll^ 
were conservative; with the scalings found in our numerical study, 
we might expect the resulting vorticity to be enhanced. This will 
form the focus of a future study. In summary, we have shown that 
the non-adiabatic pressure perturbation is non-zero in the early uni- 
verse for generic initial conditions and matter content, even if the 
intrinsic entropy perturbations are negligible. 
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